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Two  methods  are  used  to  analyze  ELF  propagation  near  lateral  nonhomogeneities 
in  the  earth-ionosphere  waveguide.  The  first,  a  full-wave  integral  equation,  is 
accurate  but  expensive  to  solve,  and  the  second,  a  two-dimensional  ray  trace,  is 
only  seraiquantitative  but  provides  physical  insight.  Both  methods  are  applied 
to  several  model  SPEs,  including  one  hased  on  the  23  November  1982  event,  which 
is  of  special  interest  because  simultaneous  measurements  of  ELF  signal  anomalies 
and  incident  proton  fluxes  are  available.  Calculations  predict  a  region  of 
weak  field  near  the  Gulf  of  Alaska,  where  a  submarine- borne  receiver  measured  an 
unusually  severe  signal  loss.  That  behavior  is  caused  by  lateral  refraction, 
which  bends  the  signal  away  from  the  gulf  and  into  the  disturbed  polar  cap,  where 
the  phase  velocity  of  the  TEM  mode  is  lowest.  The  theory  also  predicts  a  region  of 
strong  field  lust  inside  the  polar  cap  boudarv,  but  no  data  are  available  to  test  that 
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19.  ABSTRACT  (Continued) 

result.  Localized  regions  of  weak  or  strong  fields  will  occur  on  links 
where  the  great-circle  path  is  nearly  tangent  to  a  boundary  between 
large  zones  with  different  ionospheric  or  ground  properties  and  hence, 
different  waveguide  phase  velocities. 


SUMMARY 


-V 

Two  methods  are  used  to  analyze  extremely  low  frequency  (ELF) 
propagation  near  lateral  nonhomogeneities  in  the  earth-ionosphere 
waveguide.  The  first,  a  full-wave  integral  equation,  is  accurate 
although  expensive  to  solve  and  the  second,  a  two-dimensional  ray 
trace,  is  only  semiquantitative  but  provides  physical  insight.  Both 
methods  are  applied  to  several  model  solar  proton  events  (SPEs), 
including  one  based  on  the  23  November  1982  event,  which  is  of  special 
interest  because  simultaneous  measurements  of  ELF  signal  anomalies  and 
incident  proton  fluxes  are  available.  Calculations  predict  a  region  of 
weak  field  near  the  Gulf  of  Alaska,  where  a  submarine-borne  receiver 
measured  an  unusually  severe  signal  loss.  That  behavior  i3  caused  by 
lateral  refraction,  which  bends  the  signal  away  from  the  gulf  and  into 
the  disturbed  polar  cap,  where  the  phase  velocity  of  the  TEM  mode  is 
lowest.  The  theory  also  predicts  a  region  of  strong  field  just  inside 
the  polar  cap  boundary,  but  no  data  are  available  to  test  that  result. 
Localized  regions  of  weak  or  strong  fields  will  occur  on  links  where 
the  great-circle  path  is  nearly  tangent  to  a  boundary  between  large 
zones  with  different  ionospheric  or  ground  properties  and  hence, 
different  waveguide  phase  velocities. 
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I.  INTRODUCTION 


Extremely  low  frequency  (ELF)  signals  radiated  from  the  Wisconsin 
Test  Facility  (WTF)  often  exhibit  anomalies  that  are  too  strong  and 
localized  to  be  caused  by  global  changes  in  the  attenuation  rate, 
which  is  low.  Moreover,  those  anomalies  could  not  be  caused  by  mode 
interference,  because,  at  ELF,  only  the  TEM  mode  can  propagate.  A 
satisfactory  explanation  requires  a  theory  that  accommodates  lateral 
inhomogeneities  in  the  earth-ionosphere  waveguide. 

It  is  not  surprising  that  waveguide  nonstratification  must  be 
accounted  for  at  ELF.  Even  such  large  inhomogeneities  as  sporadic-E 
patches,  the  disturbed  polar  cap,  and  the  day/ night  terminator  can 
cause  the  properties  of  the  earth-ionosphere  waveguide  to  change 
markedly  over  the  huge  wavelength  or  Fresnel  zone  of  an  ELF  signal.  At 
higher  frequencies,  the  wavelength  or  Fresnel  zone  is  usually  much 
smaller  than  the  scale  lengths  of  lateral  variations  in  the  waveguide. 

Field  and  Joiner  [1979,  1982]  derived  an  integral  equation  that 
describes  ELF  fields  when  the  earth- ionosphere  waveguide  is  not 
stratified.  Since  no  directly  applicable  propagation  data  were  avail¬ 
able  at  that  time,  they  applied  the  theory  to  hypothetical  ionospheric 
disturbances,  to  nominal  models  of  the  day/night  terminator,  and  to 
the  disturbed  polar  cap.  Their  results  show  that  strong  localized 
anomalies  could  be  caused  by  focussing,  diffraction,  and  reflection  of 
the  TEM  mode. 

Pappert  [1980]  used  an  integral  formulation  to  analyze  the  effects 
of  sporadic-E  patches  on  ELF  propagation.  It  was  found  that  at  night, 
such  patches  could  cause  severe  fades.  The  physical  mechanism  is 
resonant  attenuation  that  occurs  when  the  vertical  wavelength  of  the 
ELF  wave  matches  the  optical  thickness  of  the  sporadic-E  layer.  Be¬ 
cause  of  the  masking  effect  of  the  D-region,  that  resonant  attenuation 
would  not  be  seen  in  the  daytime. 

Bannister  [1982]  summarizes  nocturnal  ELF  anomalies  measured 
simultaneously  in  the  northeastern  U.S.  and  on  board  submarines  in  the 
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North  Atlantic.  The “TOrtheas tern  U.S.  signals  faded  by  many  decibels 
in  a  few  hours  although  the  propagation  paths  were  only  about  1.6  Nfin 
long.  In  the  North  Atlantic,  three  submarines  separated  by  less  than 
2  Mm  measured  signals  that  differed  from  one  another  by  up  to  7  dB. 
Some  of  those  anomalies  had  forms  similar  to  ones  measured  in  the 
northeastern  U.S.  a  few  hours  earlier.  Moreover,  they  exhibited 
amplitude  and  phase  fluctuations  consistent  with  those  predicted  by 
Pappert  Cl 980].  Therefore,  they  were  likely  caused  by  nocturnal 
sporadic-E  patches  drifting  from  west  to  east.  A  number  of  anomalies 
were  observed  in  the  northeastern  U.S.,  but  no  simultaneous  measure¬ 
ments  were  made  at  other  sites.  Although  sporadic-E  patches  could 
have  caused  those  anomalies,  an  alternative  explanation  is  standing- 
wave  or  diffraction  patterns  caused  by  the  interaction  of  the  TEM  mode 
with  the  polar  cap  boundary. 

Katan  and  Bannister  [forthcoming]  also  report  several  anomalies 
measured  by  the  Naval  Underwater  Systems  Center  (NUSC)  at  a  number  of 
locations  during  the  solar  proton  events  (SPEs)  of  13  February  1978, 

22  to  26  November  1982,  and  8  December  1982.  The  measurements  of  the 
WTF  signal  made  in  Connecticut  and  on  board  submarines  in  the  Gulf  of 
Alaska  and  the  North  Atlantic  are  the  most  relevant  to  this  report.  In 
those  cases  anomalies  were  observed,  although  the  great-circle 
propagation  path  passed  near,  but  not  beneath,  the  disturbed  polar 
cap.  During  the  weak  23  November  1982  SPE  a  submarine  in  the  Gulf  of - 
Alaska  measured  a  fade  deeper  than  that  occurring  on  longer  paths  that 
passed  through  the  main  portion  of  the  disturbed  polar  cap  during  much 
stronger  events.  That  behavior  cannot  be  explained  with  the  often  used 
WKB  treatment  of  ELF  propagation  which  attributes  all  propagation 
phenomena  to  the  state  of  the  ionosphere  directly  over  the  great- 
circle  path. 

This  report  treats  as  a  test  case,  the  anomalies  measured  in  the 
Gulf  of  Alaska  during  the  22  to  26  November  1982  SPE.  It  compares  the 
results  obtained  using  two  computational  methods:  (1)  a  full-wave 
integral  equation  solution,  and  (2)  a  two-dimensional  ray  trace  that 
accounts  for  lateral  refraction  of  the  TEM  mode  in  the  nonstratif ied 
earth- ionosphere  waveguide. 


The  N0AA6  and  G0ES2  satellites  measured  incident  proton  fluxes  as 
a  function  of  time  and  latitude  during  the  above  three  SPEs.  Those 
important  auxiliary  data  are  used  as  inputs  to  air-chemistry  codes  to 
calculate  the  electron  and  ion  density  height-profiles  at  the  time  of 
the  ELF  anomalies.  Those  profiles  are,  in  turn,  used  as  inputs  to  our 
propagation  calculations.  Thus,  the  measured  and  calculated  ELF  fields 
can  be  compared  under  disturbed  conditions  using  actual,  rather  than 
nominal,  ionospheric  profiles. 


II.  ELF  ANOMALIES  MEASURED  DURING  SOLAR  PROTON  EVENTS 


Measurements  of  ELF  amplitude  and  phase  were  made  by  Katan  and 
Bannister  [forthcoming]  during  several  SPEs  that  occurred  between  the 
years  1976  and  1982.  Some  of  the  measurements  were  made  on  land,  and 
others  were  made  on  board  submarines  whose  locations  cannot  be 
specified  precisely.  However,  the  receiver  locations  are  known  well 
enough  for  the  present  analysis. 

The  strengths  of  the  SPEs  during  the  time  ELF  measurements  were 
made  varied  widely — the  strongest  on  13  February  1978  caused  8  dB  of 
riometer  absorption  and  the  weakest  on  23  November  1982  caused  only 
0.8  dB  riometer  absorption.  Care  must  be  used  in  interpreting  those 
strengths  because  riometer  absorption  and  ELF  propagation  are  governed 
by  different  regions  of  the  ionosphere.  It  is  not  unusual  for  larger 
ELF  effects  to  be  observed  during  a  weak  SPE  than  during  a  strong  one. 

This  report  concentrates  on  a  weak  SPE  that  occurred  on  23  Novem¬ 
ber  1982,  and  caused  propagation  effects  greater  than  can  be  explained 
in  terms  of  propagation  in  a  nearly  stratified  earth-ionosphere 
waveguide.  On  23  November  1982  the  WTF  transmitted  a  signal  that  was 
received  on  board  a  submarine  in  the  Gulf  of  Alaska.  The  approximate 
geometry  of  the  great-circle  propagation  path,  the  disturbed  polar 
cap,  and  the  first  Fresnel  zone  at  a  frequency  of  76  Hz  are  given  in 
Fig.  1 .  Note  that  the  great-circle  propagation  path  is  nearly  tangent 
to  the  edge  of  the  polar  cap,  which  covers  about  one-half  the  Fresnel 
zone. 

The  signal  received  in  the  Gulf  of  Alaska  on  23  November  1982  is 
plotted  in  Fig.  2.  Also  plotted  in  Fig.  2  is  an  average  ambient  signal 
for  that  location.  The  amplitude  of  the  disturbed  signal  was  3  to  4  dB 
below  the  ambient  signal  during  the  night,  and  about  2  dB  below  the 
ambient  signal  during  the  day.  This  behavior  occurred  because  the 
incident  proton  flux  had  diminished  considerably  by  sunrise  on  23 
November.  The  relative  phase  of  the  signal  remained  close  to  its 
ambient  value  throughout  the  measurement  period. 


Fig.  1— Geometry  of  propagation  path  from  WTF  to  Gulf  of  Alaska 
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-Signal  received  in  Gulf  of  Alaska: 
23  November  1982 


The  large  magnitude  of  the  SPE-induced  signal  loss  revealed  in 
Fig.  2  is  surprising  given  the  following  three  factors:  (1)  the  weak¬ 
ness  of  the  23  November  SPE,  (2)  the  relative  shortness  of  the 
propagation  path  (3*5  Mm),  and  (3)  the  great-circle  path  having  missed 
the  main  portion  of  the  disturbed  polar  cap.  For  comparison,  we  note 
that  smaller  signal  losses  were  caused  by  the  stronger  8  December  1982 
SPE  (2  dB  rlometer  absorption)  on  longer  paths  that  traversed  the 
center  of  the  polar  cap  [Katan  and  Bannister,  forthcoming].  The  ex¬ 
planation  lies  in  the  propagation  geometry  given  in  Fig.  1 ,  which 
causes  lateral  refraction  and  reflection  and,  thus,  shadow  zones  and 
interference  patterns  [Field  and  Joiner,  1982]. 
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III.  ELECTRON  AND  ION  DENSITY  PROFILES  UNDER  POLAR  CAP 

The  NOAA6  low-altitude  polar-orbiting  satellite  measured  integral 
proton  fluxes  in  four  channels  defined  by  E  >  2.5,  16,  36,  and  80  Mev, 
where  E  is  the  proton  energy.  Those  data  are  shown  in  Fig.  3,  for  the 
end  of  1982  November  [Sauer,  1983].  The  event  began  at  approximately 
1600  UT  on  23  November  and  lasted  at  least  until  the  middle  of  24 
November.  A  second  and  larger  SPE  commenced  on  26  November.  The 
fluxes  plotted  in  Fig.  3  are  averaged  over  magnetic  latitudes  higher 
than  70  deg.  The  proton  flux  deviates  only  slightly  from  its  spatial 
average,  provided  the  latitude  exceeds  60  deg.  Therefore,  the  fluxes 
given  in  Fig.  3  represent  those  throughout  the  main  portion  of  the 
disturbed  polar  cap. 

The  proton  fluxes  are  used  to  calculate  the  energy  deposition  and 
the  profiles  of  ion-pair  production  rates  in  the  ionosphere  during  the 
SPE.  These  production  rates  are  inserted  Into  the  lumped-parameter 
ionization-balance  equations,  which  are  solved  numerically  for  the 
electron  and  ion  density  profiles.  This  calculation  of  particle  den¬ 
sities  is  described  in  Appendix  A. 

The  calculated  electron  and  ion  densities  at  0800  UT  on  23  Novem¬ 
ber  (which  is  local  night  in  the  Gulf  of  Alaska)  are  presented  in 
Fig.  4.  The  profiles  calculated  between  0200  and  1100  UT  differed  only 
slightly  from  those  shown  in  Fig.  4  and  are  not  included  here.  In 
Fig.  5,  profiles  calculated  at  1800  UT  (which  is  local  daytime  in  the 
Gulf  of  Alaska)  are  presented.  Only  small  variations  were  found  among 
profiles  calculated  at  different  daylight  hours.  For  comparison. 

Figs.  4  and  5  plot  nominal  ambient  profiles. 


Proton  flux  (protons/cm 


Fig.  3— Integral  proton  fluxes  averaged  over  geo 
magnetic  latitudes  greater  than  70  deg 
(N0A66  satellite  measurement) 


W*"  •*'*>*’■  -*  »"*■"*  *  **  -*"•  * '  ■  -  "  *'*  ■  *'  .*■  .  ■■  *"*  '■  ■  ■  *'■  . 


IV.  PROPAGATION  MODEL  OF  NOVEMBER  1982  SPE 


« 


This  section  demonstrates  that  the  ELF  anomalies  shown  in  Fig.  2 
are  too  large  to  be  explained  by  increased  attenuation  along  the 
great-circle  propagation  path.  It  then  uses  measured  proton  fluxes  to  * 

develop  a  model  of  the  nonstratif led  polar  cap  that  is  employed  in 
Secs.  V  and  VI. 

Propagation  in  Stratified  Waveguide  * 

The  equations  governing  ELF  propagation  for  laterally  uniform 
conditions  are  available  from  many  sources,  e.g.  Galejs  [1972].  To 
define  the  notation  and  illustrate  the  key  dependences,  the  equation 
for  the  spatial  dependence  of  the  fields  is  recapitulated.  The  expres-  ( 

sion  for  the  horizontal  magnetic  intensity  H  is 

-  A  (AtAr)1/2($  d)"1/2  exp(-£fi-£  d)  e‘Bd/8*7  cos*  A/m,  (1) 

t 

where  A  depends  on  the  antenna  moment,  frequency,  and  ground  conduc¬ 
tivity,  but  not  on  the  ionosphere;  X  is  the  wave  length  in  megameters; 
and  <p  is  the  angle  between  the  propagation  path  and  the  end-fire 
direction.  The  excitation  factors  are  AT  and  AR  at  the  transmitter  and  1 

receiver  locations,  c/v  is  the  relative  phase  velocity,  3  is  the 
attenuation  rate  in  decibels  per  megaraeter  of  propagation,  and  d  is 
the  distance  from  the  source  in  megameters. 

The  amplitude  of  a  long-range  ELF  signal  is  most  sensitive  to  1 

changes— first  in  the  attenuation  rate  S,  and  second,  in  the  excita¬ 
tion  factor  A.  The  phase  depends  mainly  on  c/v. 

In  the  idealized  limit  of  a  perfectly  reflecting,  sharply  bounded 
ionosphere  at  a  height -Hq  above  the  ground,  the  excitation  factor  A  is 
equal  to  1/Hq.  Even  for  diffuse  ionospheric  boundaries,  as  treated  in 
this  report,  the  magnitude  of  A  is  of  the  same  order  as  the  reciprocal 
of  the  nominal  ionosphere  reflection  heights.  Thus,  an  ionospheric 
disturbance  that  lowers  the  effective  height  will  increase  A. 


The  ELF  data  shown  in  Fig.  2  compare  the  amplitude  and  phase  of 
the  signal  during  the  23  November  SPE  with  the  amplitude  and  phase  of 
the  signal  under  average  undisturbed  conditions.  Therefore,  to  avoid 
complexity,  the  difference  between  the  ambient  and  disturbed  signals 
was  computed  instead  of  the  signals  themselves. 

Equation  (1)  can  be  used  to  derive  the  simple  expression  given 
below  for  the  SPE-induced  amplitude  changes  AH^,  and  phase  changes  A0. 
This  can  be  done  if  we  ignore  minor  phase  variations  in  A  and  c/v 
which,  strictly  speaking,  have  complex  values  with  small  imaginary 
parts,  and  if  we  ignore  the  minor  dependence  of  amplitude  on  c/v. 

AH^  •  (8  -  8SPE)D  +  10  log  (a|PEa|PE/AtAr)  dB,  (2) 

A0  -  91 .7 (c/v  -  (c/v)SPE]  D  deg.  (3) 

Quantities  corresponding  to  SPE-disturbed  conditions  are  indicated  by 
the  superscript  SPE.  The  length  of  the  propagation  path  that  lies 
beneath  the  disturbed  cap  is  given  by  D.  Equations  (2)  and  (3)  apply 
when  the  earth-ionosphere  waveguide  behaves  as  though  it  were  nearly 
stratified. 

Calculated  Propagation  Parameters 

The  literature  supplies  full-wave  methods  for  calculating  B,  c/v, 
and  A  for  virtually  any  ionospheric  height  profile,  as  well  as  numeri¬ 
cal  results  for  many  models  of  ambient  and  disturbed  ionospheres 
[Budden,  1961a;  Field,  1970;  Wait,  1970;  Galejs,  1972;  Pappert  and 
Moler,  197^;  Greifinger  and  Greifinger,  1978].  The  profiles  of  Figs.  H 
and  5  were  used  as  input  to  the  Pacif ic-Sierra  Research  Corporation 
(PSR)  full-wave  computer  code  [Field,  1970].  This  code  calculates  the 
propagation  parameters  under  the  center  of  the  polar  cap  and  under 
ambient  conditions.  Table  1  lists  the  results  of  those  calculations. 

Equations  (2)  and  (3)  assume  a  nearly  stratified  waveguide  and  do 
not  apply  to  the  geometry  shown  in  Fig.  1.  However,  they  can  be  used 
to  obtain  an  upper  bound  on  the  anomalous  attenuation  that  the  signal 
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Table  1 


Calculated  propagation  parameters  at  76  Hz 


Date 

Time  (UT) 

8(dB/.Mm) 

c/v 

Atm"1) 

Ambient 

Night 

0.87 

1.14 

7.5  x  io"6 

23  November 

1982 

0800  (night) 

1.62 

1.24 

9.2  x  10"6 

Ambient 

Daylight 

0.94 

1.19 

9.2  x  10"6 

23  November 

1982 

1600  (daylight) 

1.28 

1.24 

1.0  x  io"5 

could  have  suffered  propagating  on  a  great-circle  path  from  the  WTF  to 
the  Gulf  of  Alaska.  To  obtain  that  bound,  only  the  first  term  on  the 
right  side  of  Eq.  (3)  is  retained.  This  is  equivalent  to  the  assump¬ 
tion  that  the  entire  path  is  beneath  the  disturbed  polar  cap  while  the 
terminals  arc  under  ambient  conditions.  We  overstate  the  exposed  path 
length  and  ignore  possible  enhancement.  This  enhancement  is  caused  by 
the  ionosphere  being  lowered  over  either  terminal,  thus  increasing  the 
excitation  factors. 

Insertion  of  the  attenuation  rates  from  Table  1  and  the  value  D  - 
3.5  Mm  into  Eq.  (2)  gives  AH^  -  -2.6  dB  and  AH^  -  -1.2  dB  under  night 
and  daylight  conditions,  respectively.  These  calculated  signal  losses, 
which  are  upper  bounds,  are  smaller  than  the  measured  values  of  3  to  4 
dB  and  2  dB  (see  Fig.  2).  In  fact,  the  discrepancy  between  the 
measured  losses  and  the  ones  calculated  from  stratif ied-waveguide 
theory  is  greater  than  indicated  by  the  figures  because  (1)  the  entire 
great-circle  path  was  not  exposed  to  the  central  polar  cap,  and  (2) 
the  terminals  were  so  close  to  the  polar  cap  boundary  that  some  in¬ 
crease  in  excitation  factor  (say,  one-half  of  the  values  shown  in 
Table  1)  would  be  expected.  If  those  adjustments  are  made,  Eq.  (2) 
predicts  virtually  no  signal  loss  in  the  daylight  and  about  1  dB  loss 
at  night.  Moreover,  Eq.  (3)  predicts  a  -32  deg  phase  shift  at  night 
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and  a  -16  dB  phase  shift  in  daylight,  but  virtually  no  phase  shifts 
were  observed. 

These  results  show  that  the  measurements  cannot  be  explained  by  a 
theory  that  omits  lateral  refraction,  reflection,  and  diffraction  at 
the  boundary  of  the  disturbed  polar  cap. 

Polar  Cap  Boundary  Model 

In  order  to  model  the  polar  cap  boundary,  measurements  of  the 
proton  flux  versus  geomagnetic  latitude  are  needed.  Because  the  data 
for  the  23  November  SPE  reveal  only  the  time  dependence  of  the  flux 
averaged  over  the  cap,  the  latitudinal  dependence  must  be  inferred 
from  data  measured  during  other  SPEs.  For  example,  Fig.  6  shows  the 
flux  versus  latitude  for  the  8  December  1982  event.  The  fluxes  are 
approximately  constant  above  60  deg  and  fall  off  at  a  fairly  constant 
rate  between  about  60  and  50  deg.  Therefore,  we  assume  the  disturbed 
polar  cap  to  be  uniform  between  the  north  geomagnetic  pole  and  60  deg 
latitude,  which  is  about  3-0  Mm  from  the  pole.  The  diffuse  boundary, 
or  transition  zone,  extends  to  50  or  55  deg  and  is  0.5  to  1  Mm  wide. 
Both  WTF  and  the  Gulf  of  Alaska  are  at  about  55  deg  north  geomagnetic 
latitude  and  are  in  the  boundary  region. 

We  assume  a  flat  earth,  an  isotropic  ionosphere,  symmetry  about  the 
north  geomagnetic  pole,  and  model  the  propagation  constant  as: 

S(x ,  y)  -  SSPE  +  (SamB  “  SspE)/{1  +  exp[-(r  -  rQ)/5r]},  (4) 

where  S  is  related  to  the  phase  velocity  and  attenuation  rate  by  the 
formulas  c/v  »  reS,  and  8  ■  -8.6k  imS.  In  addition,  SAMB  and  SspE 
denote  the  values  of  S  under  ambient  and  disturbed  conditions,  respec¬ 
tively.  The  distance  from  the  pole  is  r,  rg  is  the  radius  of  the 
center  of  the  polar  cap  boundary,  and  5r  is  a  scale  distance  that  we 
use  below  to  define  the  width  of  the  boundary.  Figure  7  diagrams  this 
model  and  defines  some  parameters. 

If  the  boundary  width  Ar  is  defined  as  the  distance  over  which  S 
undergoes  95  percent  of  its  transition  from  disturbed  to  ambient,  and 
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Fig.  7— Model  of  polar  cap  propagation  constant 

r1  and  p2  are  the  inner  and  outer  radii  of  the  boundary  region, 
respectively,  it  follows  that  (see  Appendix  B) 

r2  -  r1  ♦  Ar  ,  (5) 

r0  “  (r1  +  r2)/2  • 

Ar  -  7.3  5r  .  (7) 

The  transmitter-to-receiver  path  length  is  3.5  Mm  on  a  great-circle 
path.  We  will  assume  that  the  Gulf  of  Alaska  receiver  is  approximately 
the  same  distance  from  the  north  geomagnetic  pole  as  the  WTF  transmit¬ 
ter.  Since  the  exact  location  of  the  receiver  is  not  certain,  we  will 
regard  the  gulf  as  a  1  Mm  square  area,  the  center  of  which  is  3.5  Mm, 
from  both  the  WTF  and  the  north  geomagnetic  pole.  We  will  assign  r-j 
the  value  of  3.0  Mm, which  allows  the  model  of  the  propagation  medium 
in  Fig.  8. 


V.  TWO-DIMENSIONAL  RAY  TRACE:  LATERAL  REFRACTION  OF  TEM  MODE 


The  WTF  and  Gulf  of  Alaska  are  at  about  55  deg  geomagnetic  north 
latitude  and  probably  lie  in  the  boundary  region,  where  strong 
transverse  gradients  in  S  will  refract  the  field.  More  simply,  the 
propagation  constant  S  behaves  as  a  refractive  index  and  tends  to  bend 
the  fields  toward  the  inner  cap,  where  the  phase  velocity  is  slowest. 
Under  the  proper  conditions,  this  effect  could  create  a  shadow  zone  to 
which  rays  emanating  from  the  WTF  cannot  gain  access. 

To  test  this  hypothesis  a  simple  ray  trace  is  used.  It  is  sum¬ 
marized  in  this  section  and  derived  in  Appendix  B.  The  results  are 
only  semi quantitative ,  because  the  validity  criterion  for  ray  tracing 
is  marginally  satisfied  for  most  of  the  assumed  boundaries  and  is 
actually  violated  for  the  abrupt  ones.  Therefore,  ray  tracing  wilj. 
overstate  the  refraction.  In  Sec.  VI,  the  fields  are  recalculated 
using  a  full-wave  integral  equation. 

We  assume  a  two-dimensional  propagation  medium,  whose  propagation 
constant  S  is  given  by  Eq.  (4).  Because  the  propagation  constant  is 
analogous  to  a  refractive  index,  the  problem  is  mathematically  identi¬ 
cal  to  tracing  rays  that  are  obliquely  incident  on  an  isotropic  iono¬ 
sphere  that  varies  in  two  dimensions,  but  is  uniform  in  the  direction 
perpendicular  to  the  plane  of  propagation.  To  isolate  refractive 
effects,  we  ignore  the  imaginary  part  of  S  and,  hence,  absorption. 
Signal  anomalies  calculated  in  this  fashion  are  caused  solely  by 
lateral  focusing  or  defocus ing  of  the  TEM  mode. 

The  following  geometric  parameters  are  assumed: 


r1  -  3  Mm  , 

(8) 

Ar  ■  1  Mm  . 

(9) 

Rays  are  traced  for  a  total  of  six  levels  of  disturbance,  including: 
(1)  weak,  equivalent  to  the  23  November  1982  SPE,  characterized  by 
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re^SPE  ”  1*25,  (2)  moderate,  a  bit  stronger  than  the  8  December  1982 
SPE,  characterized  by  reSSPE  -  1.35,  and  (3)  very  strong,  charac¬ 
terized  by  reSgpg  -  1.50.  In  each  case  ambient  conditions  are  charac¬ 
terized  by  reSAM B  -  1.15,  which  corresponds  to  nighttime  propagation. 

A  top  view  of  the  ray  trajectories  is  diagrammed  in  Fig.  9.  In  all 
cases  the  rays  bend  toward  the  polar  cap  and  become  less  dense  in  the 
southern  Gulf  of  Alaska,  shown  as  the  shaded  area.  The  signal  is 
weakest  where  the  rays  are  least  dense  and  strongest  where  they  are 
densest.  As  is  to  be  expected,  the  refraction  increases  as  the  SPE 
becomes  stronger.  For  all  but  the  weakest  of  the  modeled  events,  a 
caustic  is  formed  near  the  inner  edge  of  the  boundary  region. 

Figure  10  reveals  the  dependence  of  the  ray  trajectories  on  the 
thickness  of  the  boundary.  Those  results  were  calculated  for  the  weak 
SPE  (S  -  1.25)  and  boundary  thicknesses  ranging  from  0.5  to  1.0  Mm.  As 
before,  the  radius  of  the  inner  polar  cap  is  assumed  to  be  3.0  Mm.  The 
results  given  in  Figs.  9  and  10  show  that  the  signal  depends  strongly 
on  the  boundary  thickness  as  well  as  on  the  strength  of  the  distur¬ 
bance. 

The  effects  of  an  SPE  on  3ignal  pattern  (shown  in  Figs.  9  and  10) 
would  be  azimuthally  symmetric  if  the  ionosphere  were  stratified. 
Figures  9  and  10  indicate  the  difference  between  the  actual  field  in  a 
region,  and  what  would  have  been  received  if  the  WTF  radiation  pattern 
were  not  distorted  by  lateral  refraction.  Although  the  ray  trace  is  a 
convenient  means  of  estimating  the  locations  and  relative  strengths  of 
signal  concentrations  and  rarefactions.  Figs.  9  and  10  cannot  be  used 
to  compare  ambient  fields  with  those  measured  during  an  SPE.  That 
calculation  would  require  each  ray  to  be  weighted  by  a  number  that 
accounts  for  the  waveguide  excitation  factor. 

Figures  9  and  10  show  that  lateral  refraction  can  be  substantial. 

£ 

It  is  estimated  that,  for  a  1-Nta  thick  boundary,  refraction  changes 
the  signal  in  the  southern  Gulf  of  Alaska  -8  dB  relative  to  a  sym¬ 
metric  signal  for  the  strong  SPE;  -3.5  dB  for  the  moderate  SPE;  and 

ft 

The  received  field  i3  inversely  proportional  to  the  square  root  of 

the  distance  between  rays. 


-1.4  dB  for  the  weak  SPE.  These  effects  increase  as  the  boundary 
becomes  more  abrupt.  For  the  weak  SPE,  which  corresponds  to  the  23 
November  event,  the  ray  trace  predicts  that  defocusing  in  the  southern 
Gulf  of  Alaska  is  as  much  as  -6  dB  for  the  narrowest  boundary  (0.5  Mm) 
and  -2  dB  for  a  boundary  width  of  0.7  Mm.  Note  that  these  refractive 
effects  are  over  and  above  losses  attributable  to  attenuation. 

The  field  depends  strongly  on  receiver  location.  For  all  cases 
shown,  a  region  of  intensification  (focusing)  occurs  just  a  few 
hundred  kilometers  north  of  the  region  of  minimum  field.  Since  the 
exact  location  of  the  submarine-borne  receiver  is  unknown,  a  detailed 
comparison  cannot  be  made  between  experiment  and  theory.  The  concept 
of  energy  refracting  away  from  the  Gulf  of  Alaska  and  into  the  polar 
cap  appears  consistent  with  the  measured  SPE-induced  signal  loss  that 
was  too  severe  to  have  been  caused  by  attenuation  alone. 


Fig.  10— Ray  tracing  for  weak  SPE  and  various  transition  zone  widths 


VI.  FULL-WAVE  CALCULATIONS  USING  INTEGRAL  WAVE  EQUATION 


This  section  presents  results  calculated  from  an  integral  equation 
developed  by  Field  and  Joiner  Cl  979]  to  analyze  ELF  propagation  under 
conditions  where  the  earth-ionosphere  waveguide  is  not  stratified. 

That  full-wave  equation  accounts  for  a  number  of  phenomena  that  the 
ray  trace  omits,  including  attenuation,  diffraction,  and  lateral 
reflection.  A  major  disadvantage  of  the  equation  is  the  large  computer 
expense  needed  to  achieve  a  solution.  Therefore,  only  a  few  sample 
cases  are  shown  below.  Concentration  is  on  effects  caused  by  diffrac¬ 
tion  and  reflection. 

Integral  Wave  Equation 

Field  and  Joiner  [1979]  write  the  field  at  the  ground  as 

E  -  A  A(x,  y)  f(x,  y)  F(z)  ,  (10) 

where  A  is  a  constant  involving  dipole  moment,  wave  frequency,  and 
ground  conductivity;  A  is  the  excitation  factor  and  a  function  of 
position;  F  is  the  vertical  dependence;  and  V  is  the  lateral  depend¬ 
ence  of  the  signal.  The  lateral  dependence  can  be  found  by  solving 

[V2  +  k2S2(x,  y)]  V  -  0  ,  (11) 

where 


Here  S  is  the  propagation  constant  and  is  found  by  imposing  boundary 
conditions  on  F(z)  at  the  ground  and  in  the  ionosphere. 

The  WTF  consists  of  two  nearly  perpendicular  horizontal  dipoles. 
Katan  and  Bannister  [forthcoming]  report  that  the  phase  difference 
between  the  two  antenna  elements  was  290  deg  during  the  23  November 
tests.  That  phase  angle  gives  a  radiation  pattern  that,  for  the  pur- 
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poses  of  the  present  report,  can  be  assumed  symmetric.  If  we  further 
assume  the  transmitter  to  be  at  the  origin  of  a  cylindrical  coordinate 
system,  then  for  a  stratified  waveguide  where  S  (x,  y)  -  Sq,  the 
solution  of  Eq.  (11)  is 


H^2) (k  Sc 


r) 


(12) 


where  r  is  the  distance  from  the  origin,  and  H  is  the  Hankel  function. 

Strictly  speaking,  the  separation  of  variables  that  leads  to  the 
solution  given  by  Eq.  (12)  is  invalid  when  the  ionosphere  is  not 
stratified.  However,  under  most  conditions  the  eigenvalue  S  is 
governed  primarily  by  the  local  ionosphere.  That  assumption  allows  the 
solution  to  Eq.  (10)  to  be  written  in  the  form 

¥  -  W(x,  y)  H£2)(k  SAMB  r)  ,  (13) 


where  W(x,  y)  is  a  propagation  factor  that  is  unity  under  undisturbed 
conditions  and  gives  the  relative  effect  of  a  nonstratif led  distur¬ 
bance  on  the  signal.  A  lengthy  derivation  [Field  and  Joiner,  1979] 
reveals  that  W(x,  y)  satisfies  the  following  integral  wave  equation: 


W(x,  y) 


1 


dy' 


s  (x,  y) 


,(2) 


k  S 


i)  h°2)( 


AMB 


H<2)(k 


k  S 


AMB  lJ 


r) 


W(x\  y') 


(14) 


where 


2  2  2 
r  «  x  +  y 


r2  -  (x')2  +  (y')2  , 


r|  -  (x  -  x')2  +  (y  -  y')2 
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Numerical  Results 

As  was  done  for  the  ray-trace  calculations  presented  in  Sec.  IV, 
the  propagation  constant  S  is  given  by  Eq.  (4),  and  only  reS  is 
retained  which  is  tantamount  to  neglecting  attenuation.  Any  SPE- 
induced  signal  changes  calculated  in  this  manner  are  caused  by  refrac¬ 
tion,  reflection,  or  diffraction  at  the  boundary.  The  changes  are  over 
and  above  losses  caused  by  attenuation.  The  results  given  below  are 
more  accurate  than  those  given  in  the  previous  section,  because  the 
ray  trace  is  not  fully  valid  for  all  boundary  thicknesses  treated,  and 
it  breaks  down  near  caustics. 

Contour  plots  of  the  change  in  W  caused  by  a  weak  SPE  having 
boundary  thicknesses  of  0.5  and  1  Mm,  respectively,  are  illustrated  in 
Figs.  11  and  12.  Those  parameters  represent  the  23  November  event.  The 
figures  show,  in  decibels,  the  signal  during  the  event  relative  to  the 
signal  during  ambient  night.  Figure  13  is  analogous  to  Figs.  11  and 
12,  but  applies  to  a  strong  SPE  with  a  1-Mm  boundary  thickness.  Note 
that  for  a  stratified  ionosphere,  the  contours  would  be  azimuthally 
symmetric  about  the  transmitter. 

The  results  are  qualitatively  similar  to  the  ray  traces.  In  all 
cases  there  is  a  signal  increase  in  the  northern  gulf  where  the  ray 
trace  produced  a  caustic,  and  a  signal  decrease  in  the  southern  gulf 
where  the  ray  trace  showed  a  decrease  in  ray  density.  Presumably,  that 
decrease  corresponds  to  the  anomalous  signal  loss  measured  during 
certain  SPEs. 

The  signal  distortions  shown  in  Figs.  11  through  13  are  not  as 
strong  as  those  given  by  the  ray  trace,  which,  for  reasons  given 
above,  overstates  the  refraction. 


VII.  CONCLUSIONS 


In  order  to  interpret  anomalous  signal  losses  measured  during 
several  SPEs ,  two  methods  were  developed  for  handling  ELF  propagation 
near  lateral  nonhomogeneities  in  the  earth-ionosphere  waveguide:  (1)  a 
full-wave  integral  equation  that  is  accurate,  but  expensive  to  solve, 
and  (2)  a  two-dimensional  ray  trace  that  provides  physical  insight, 
although  it  is  only  semi  quantitative. 

Both  methods  were  applied  to  several  model  SPEs,  including  one 
based  on  the  23  November  1982  event,  which  is  of  special  interest 
because  simultaneous  measurements  of  ELF  signal  anomalies  and  incident 
proton  fluxes  are  available.  All  calculations  indicated  a  region  of 
weak  field  in  the  southern  Gulf  of  Alaska,  where  a  submarine-borne 
receiver  measured  an  unusually  severe  signal  loss.  That  behavior  is 
caused  by  lateral  refraction,  which  distorts  the  WTF  radiation  pattern 
by  bending  energy  away  from  the  gulf  and  into  the  disturbed  polar  cap, 
where  the  phase  velocity  of  the  TEM  mode  is  lowest.  Accordingly,  the 
theory  also  predicts  a  region  of  strong  field  just  inside  the  polar 
cap  boundary,  but  no  data  are  available  to  test  that  result. 

Localized  regions  of  weak  or  strong  fields  are  by  no  means 
restricted  to  situations  where  the  great-circle  propagation  path 
grazes  the  edge  of  the  polar  cap,  as  was  the  case  for  the  WTF  to  Gulf 
of  Alaska  link  considered  here.  They  will  occur  on  any  link  where  the 
great-circle  path  is  nearly  tangent  to  a  boundary  between  large  zones 
with  different  ionospheric  or  ground  properties  and,  hence,  different 
waveguide  phase  velocities.  Examples  of  such  boundaries  include  the 
terminator  and  the  Greenland  shoreline. 


Appendix  A 


CALCULATION  OF  IONIZATION  PROFILES  FROM  PROTON  FLUXES 


In  this  appendix  we  discuss  the  calculation  (using  proton  flux 
data)  of  the  ionospheric  charged  particle  densities.  First,  the  proton 
data  are  used  to  determine  rates  of  ion-pair  production  at  altitudes 
that  affect  ELF  propagation.  Second,  those  rates  are  input  to  air- 
chemistry  equations  to  obtain  height  profiles  of  electron  and  ion 
densities. 

Calculation  of  ion-pair  production  rates  requires  knowledge  of  the 
proton  flux  in  a  large  number  of  narrow  energy  bin3  between  about  1 
and  300  MeV.  To  find  the  number  of  protons  at  energy  E  or  greater  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA)  fits  the  four 
coarse  energy  bands  of  the  spacecraft  data  to  the  following  empirical 
formula  for  the  integral  proton  flux: 

J  *  J0  e'  R/R0  protons/cm2-sec-sr  ,  (A.1) 


where 


R  -  E2  +  2EE0  MeV  .  •  (A. 2) 

In  the  above  equations,  Eq  is  the  proton  rest  energy  (936  MeV),  and  JQ 
and  Rq  define  the  strength  and  energy  rigidity,  respectively,  of  the 
integral  flux. 

The  values  of  JQ  and  Rq  supplied  by  NOAA  for  23  November  1982  are 
shown  in  Fig.  A.1.  The  strength  Jq  changes  by  a  factor  of  about  2 
during  the  daylight  and  night  periods.  The  rigidity  Rq,  however, 
changes  little.  In  fact,  even  the  models  generated  using  Rq  for  dif¬ 
ferent  times  during  daylight  and  night  show  little  difference. 
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Fig.  A. 1— Strength  and  rigidity  parameters: 
23  November  1982 


In  Fig.  A. 2  [Reid,  1978],  profiles  of  ion-pair  production  rates 
per  proton  are  presented  for  isotropic  fluxes  of  monoenergetic  protons 
of  given  energies.  To  find  the  total  production  rate,  the  number  of 
protons  of  a  particular  energy  must  first  be  found.  It  is  assumed  that 
all  the  protons  in  a  range  of  energy  have  the  same  nominal  energy.  If 
and  are  the  minimum  and  maximum  energies  of  the  range, 

respectively,  then  the  total  number  of  protons  at  the  nomimal  energy 
is  given  by 


Jtotal  ”  ^^in^  “  J^Emax^ 

Table  A.1  lists  the  energy  ranges  and  the  nominal  energies  used 


in  this  report.  Logarithmic  spacing  is  employed. 


(MeV) 


Fig.  A.2--Ionization  per  proton  of  given  energy 
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Table  A.1 


Energy  bands  used  to  find  ion-pair 
production  rates 
(In  MeV) 


Energy  Range3  Nominal  Energy 


0.6  -  2.0 

1  .0 

2.0  -  6.0 

3.0 

6.0  -  20.0 

10.0 

20.0  -  60.0 

30.0 

60.0  -  200.0 

100.0 

200.0  -  600.0 

300.0 

600.0  -  2000.0 

1000.0 

aIncludes  protons  that  deposit  their 
energy  and,  hence,  produce  ionization, 
at  altitudes  between  15  and  110  km. 


Table  A.1  and  Eq.  (A.1),  as  well  as  Fig.  A.1  allow  us  to  find  the 
production  rates  q.  These  rates  are  then  inserted  into  the  three- 
species,  lumped  parameters  deionization  model  of  Knapp  and  Schwartz 
[1975].  They  give  the  following  for  the  electron  density  Ne  and  the 
positive  ion  density  N+  at  altitude  h: 


Ne(h) 


(q  +  D) 

[A  +  D  +  adN+(h)]  ’ 


(A. 3) 


N+(h>.  f 


(A. 4) 


!:•>: 
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4ft 


1^ 


* 


%r 


£ 


Sr 


c 


where 


a 


_< 

[ao.  +  Da,  +  a 
[  i  d 

iad 

[  q  (A  +  D)"| 
_Aa  +  Da 

1/2, 

J 

f 

{ 

[a  +  D  + 

*  q  (A  +  D)' 
.Aa  +  Da  . 

1/2] 

i 

f  ’ 

(A. 5) 


In  the  above  equations  A,  D,  ,  and  ad  are  functions  of  altitude  and 
pair  production,  and  are  taken  from  figures  in  Knapp  and  Schwartz 
[1975].  Their  simple  air-chemistry  model,  developed  under  Defense 
Nuclear  Agency  auspices,  has  for  years  been  the  standard  method  of 
calculating  ionospheric  ionization  produced  by  radiation  from  high- 
altitude  nuclear  bursts.  However,  more  detailed  models  that  account 
for  large  numbers  of  ionospheric  species  do  exist,  although  the  values 
for  the  deionization  coefficients,  though  widely  used,  are  not 
precise. 

The  pair  production  rates  are  smoothed  before  insertion  into 
Eqs.  (A. 3)  through  (A. 5).  The  results  are  3hown  in  Figs.  4  and  5  on 
pp.  10  and  11.  The  particle  densities,  as  well  as  the  collision  fre¬ 
quencies,  mean  ionic  mass,  and  geomagnetic  field  strength  must  be 
specified  in  order  to  define  the  electromagnetic  properties  of  the 
ionosphere.  The  assumed  electron  and  ion  collision  frequencies  are 
diagrammed  in  Fig.  A. 3.  In  addition,  as  per  Pappert  and  Moler  [1974], 
a  nominal  ion  mass  of  32  amu  and  a  geomagnetic  field' strength  of  0.5  G 
are  used.  It  is  assumed  that  the  propagation  path  is  ea3t  to  west  and 
that  the  magnetic  field  dip  is  -80  deg.  Figures  4  and  5  also  show  the 
ambient  day  and  night  profiles  [Pappert  and  Moler,  1974]. 


--Electron  and  ion  collision  frequencies 
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Appendix  B 

RAY  TRACING  IN  TWO  DIMENSIONS 


* 


In  this  appendix  we  derive  the  ray-tracing  equations  (for  the  two- 
dimension  wave)  that  solve 

[V2  +  k2S2(x,  y)l  i|>  -  0  .  (B.l) 

The  derivations  are  based  on  Budden's  [1961b]  theory.  That  theory  is 
presented  in  full  in  order  to  indicate  the  approximations  involved. 

Derivation  of  Ray-Tracing  Equations 

The  propagation  constant  S  is  written  as  follows: 

S  =■  y  +  iX  • 


# 


The  imaginary  part  in  what  follows  is  neglected,  so  Eq .  (B.l)  becomes 


(V2  +  k2y2(x,  y)l  <P  =  0  . 


(B.2) 


Thus  y  is  equivalent  to  the  refractive  index,  and  it  will  be  referred 
to  as  such  for  convenience. 

The  solution  to  Eq.  (B.l)  can  be  written  as 


<0 


-ie(x,y) 

e 


( B .  3 ) 


In  Eq.  (B.3)  e(x,  y)  is  the  eikonal  function,  which  we  assume  can  be 
written 


e(x,  y) 


(  B .  4 ) 


where  P^(x,  y)  ,  P^(x,  y)  are  assumed  to  be  proportional  to  the  direc 
tion  cosines  of  the  two-dimensional  field's  wave  normal,  and 


2  2  2 

P  +  P  =  y 
x  y 


In  order  for  ray  tracing  to  be  possible  it  must  be  assumed  that  the 

medium  is  varying  slowly  enough  so  that  the  spatial  derivatives  of 

P  and  P  are  small  enough  that  they  can  be  neglected.  (This  point 
x  y 

will  be  discussed  later.)  If  that  assumption  is  valid,  then 


1  3e  . 

7-  and 

k  dx 


1  3e 
k  3y  ’ 


(The  subscript  T  will  hereafter  be  dropped,  because  the  existence  of 
P^,  which  is  always  zero,  is  assumed  since  f.  is  not  a  function  of  z.) 
From  Eq.  (B.7)  it  follows  that: 


V  x  p  =  0  , 


From  Eqs.  (B.3)  and  (B.4)  it  can  be  shown  that  P  ,  P  are  propor 

.  x  y 

tional  to  direction  cosines  of  the  wave  normal  at  the  point  x,  v.  In 
fact , 


•  ’*  -  •  *  .  •  .  -  *  f  „  «  , 

r>  *>*>  *>  - 
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(p2  +  p2)1/2  (p2  +  p2) 

'  x  y/  \  x  y/ 


1/2  » 


are  Che  direction  cosines  of  the  wave  normal.  If  the  wave  normal 

takes  all  possible  directions,  then  the  locus  of  the  point  (P  .  P  ) 

x  y 

in  "refractive  index  space"  defines  the  refractive  index  surface  G. 
The  equation  of  this  surface  may  be  written 


G(x,  y;  Px,  Py) 


2  2' 

P  +  P  . 
x  y 


U(x,  y;  P  P  ) 

a  y 


=  i  , 


(B  .9) 


where  y  is  made  a  function  of  P  ,  P  to  indicate  that  it  may  depend 

x  y 

on  the  direction  of  the  wave  normal.  The  ray  direction  is  normal  to 
the  refractive  index  surface,  and  thus  has  direction  cosines  given  by 


A  wavefront  is  defined  as  some  point  in  the  wave  with  a  partic¬ 
ular  phase,  such  as  the  wave  crest.  A  wavefront  travels  with  the 
ray  velocity  VR.  Let  x,  y  be  the  coordinates  at  which  the  wavefront 
intersects  the  ray.  Hence,  x,  y  are  the  components  of  the  ray  veloc¬ 
ity,  and  must  therefore  be  proportional  to  3G/9P  and  9C/3P  .  Buddon 

XV 

[1961b]  shows  that  the  constant  of  proportionality  is  the  speed  of 
light  c,  so 


(B . 10) 


Since 


G(x,  y;  Px,  Py>  -  1  , 


the  result  must  be 


dG  n  ,  dG  „ 
-7—  -  0  and  —  =  0 
dx  dy 


So,  from  the  left-hand  equation  above,  the  following  is  derived 


3G  _3G_  _3G_ 

3x  3P  3x  3P  3x 

x  y 


Using  Eq.  (B.10)  yields 


ar  t  /  3P  3P  \ 

3£  +  I  x  —  +  v  — * 
3x  c  \  3x  ^  3x  / 


And  doing  the  same  with  Eq.  (B.8)  gives 


.  /  3P  3P  \ 

c  (x  3x  y  3y  / 


which  can  also  be  written  as 


3G  1  dPx 
3x  c  dt 


Therefore, 


3x  ’ 


dy  ' 


Equations  (B.10),  (B.12),  and  (B.13)  are,  taken  all  together,  t 
canonical  equations  of  the  ray. 
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The  variable  9,  which  is  the  angle  the  wave  normal  makes  with 


the  y-axis,  is  introduced  into  the  calculations  as  follows: 


P  =  u  sin  9  , 

x 


(B. 14) 


P  =  u  cos  9 

y 


(B . 15) 


Therefore,  using  Eqs.  (B.12)  through  (B.15),  it  can  be  shown  that: 


dP 

x  du  „  d9 

dt~  =  dF  sin  9  +  W  cos  9  ^  = 


3*  • 


dPy  dy  .  .  .  d9 

"dF  =  dt  cos  9  -  U  sin  9  — 


3y 


Solving  for  d9/dt  gives 


d9  c 


f  Q  3G  .  Q  3G] 

[cos  9  ^  -  sin  9  97 J 


Differentiating  Eq .  (B.9)  yields 


(Pl  +  Py)J 


3G  _  _3_  _\  _x _ 

3x  3x  u 


I  Ik 

U  3x  ’ 


2  2 

P  +  P 

l  X  V) 


I  Ik 

U  3v  ’ 


V-V-  .  •  .N  .•*  v- 
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3G 

3P 

y 


cos  9 


U 


sin  9 
U 


4 


Equations  (B.16),  (B.17),  and  (B.18),  then,  are  the  ray  tracing 
equations : 


and , 


d0 

dt 


c 

2 

U 


[cos  0  it  .  sin  e  |tj 


(B.16) 


dx  sin  9 
dt  C  u 


(B.17) 


dy  _  cos  9 
dt  C  u 


(B.18) 


Validity  Criteria 

In  order  to  establish  the  validity  of  Eq.  (B.6)  we  will  consider 
a  situation  in  which  u  is  a  function  of  y  only  (at  least  locally). 

In  such  a  situation,  cannot  depend  on  x,  so  it  follows  that 


4 


4 


4 


4 


4 


ip  =  <pQ  exp 


-  y 

[-lk(Pxx  *j  v?)] 


which  resembles  the  equation  of  a  field  at  oblique  incidence  to  the 
ionosphere. 

Just  as  nS  is  independent  of  z  in  the  case  of  propagation  in  the 
ionosphere,  it  can  be  shown  that  P^  is  independent  of  y: 


MLJdSi 
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V2  *  ,  if*  +  £* 

7  ax2  ay2 


o  .  3P  \ 

*V^)  41 


y  k  3y 


•  -k2  v2(1+ZJ  V1) 


So,  the  solution  to  Eq .  (B.2)  is  valid  if 


2  3y 


By  differentiating  Eq .  (B.5)  and  using  Eq.  (B.15)  we  obtain 


y  =  1  3y 

3y  cos  9  3y 


Therefore,  the  following  is  the  validity  criterion: 


_ I _  1M  «  j 

,2  3y  L 

ky  cos  6 


(B. 19) 


This  will  fail  for  any  value  of  or  y  if  9  is  near  90  deg. 


Typically  Ay  =  0.1,  Ay  S  1  Mm,  y  S  1.20,  and  k  =  1.57  1/Mm 
at  76  Hz.  Thus,  at  9  =  0,  if  we  say  3y/Sy  »  Ay/Av,  then, 


1  3y 


0.044 


>‘vv«V.V  Vv  v  .•  V'aVW 


Therefore,  for  an  accurate  ray  trace, 


0.044  <<  cos  0 


or. 


8  <<  87  deg  . 


To  obtain  a  more  precise  understanding,  0.1  is  taken  as  the  upper 
limit  on  the  left-hand  side  of  Eq .  (B.19).  This  yields 


0.044  < 
cos  8 


0.1 


or. 


cos  8  >  0.44 

Therefore,  rays  that  make  angles  with  the  y-axis  of  0  <  65  deg  are 
valid.  If  0.2  is  the  upper  limit,  then  rays  with  0  <  77  deg  are 
valid . 

In  the  above,  y  was  assumed  to  change  only  in  the  y-direction. 
In  fact,  the  disturbance  of  23  November  1982  is  a  function  of  the 
distance  from  the  north  pole.  However,  Eq .  (B.3)  must  be  a  solution 
at  each  point  in  space.  Therefore  if  the  y-axis  is  taken  locally  to 
be  along  the  gradient  of  y,  then  Eq .  (B.19)  should  yield  the  local 
validity  condition,  which  is  expected  to  fail  at  right  angles  to  the 
direction  of  Vy .  The  validity  condition  can  be  written  as 


1  _1_  dy 

cos  8  2  dr 

y 


«  l 


( B .  20 


In  the  above  calculations,  3y/3r  is  approximated  by  Au/Ar.  How¬ 
ever,  using  Eq.  (4)  gives  the  value  of  3y/3r  at  all  points.  Since  no 
absorption  is  assumed  in  Sec.  V,  S  can  be  substituted  for  y.  Then, 
letting  AS  =  S  -  S  ,  Eq .  (4)  can  be  written 


-4  7- 


S(r)  -  -  f 


r  (r  “  r0)' 

) 

1  -  exp 

<5r 

■  {  +  s  +  M 

r  (r  -  ro)i 

(  SAMB  +  2 

1  +  exp 

L  Sr  J 

) 

The  inner  and  outer  radii  of  the  boundary  region,  r^  and  r^,  respec¬ 
tively,  [as  shown  in  Fig.  7  on  p.  17  and  defined  by  Eqs.  (5)  and  (6) 
on  p.  17]  can  be  expressed  by  the  relations 


S<V  -  +  f  °'95  +  SAMB+¥ 


S(r2)  -  -f  0.95  +  S^+f  . 


Solving  for  r^  and  yields 


(ri  -  r0} 


=3.66  , 


(r2  '  r0} 


=  -3.66  , 


^r2  "  rl^  = 


Therefore, 
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4 


Sr 


Ar 

7.33 


(B . 21) 


The  derivative  of  S  becomes 


dS 

dr 


(B . 22) 


This  derivative  has  a  maximum  at  r  *  rrt,  so  the  maximum  is 


dS  _  AS  J_  7.33  AS  AS 

dr  4  Sr  4  Ar  Ar 


(B.23) 


The  left-hand  side  of  Eq .  (B.20)  can  now  be  calculated  for  the  follow¬ 
ing  three  situations  discussed  in  Sec.  V:  (1)  the  weak  SPE  with  a 
1  Mm  transition  zone,  (2)  the  strong  SPE  with  a  1  Mm  transition  zone, 
and  (3)  the  weak  SPE  with  a  0.5  Mm  transition  zone.  Table  B.l  gives 
the  results. 

The  values  in  Table  B.l  were  calculated  for  the  maximum  value  of 
the  derivative.  The  validity  condition  is  improved  considerably  at 
distances  nearer  to  or  farther  from  the  north  pole  than  r^.  This  is 
illustrated  by  the  plot  in  Fig.  B.l,  which  shows  the  validity  crite¬ 
rion  for  situation  3  at  80  deg  as  a  function  of  distance  from  the 
north  pole. 


i 


Table  B.l 


9 

Validity  Criterion 

For  Three  Examples 

9 

Example  1 

Example  2 

Example  3 

Angle 

with  Respect 
to  Gradient 
(deg) 

SSPE  =  1.25 

AS  a  0.1 

Ar  a  1  Mm 

SSPE  =  1.45 

AS  a  0.3 

Ar  a  1  Mm 

SSPE  =  1.25 

AS  a  0 . 1 

Ar  S  0.5  Mm 

0 

0 

0.081 

0.207 

0.162 

60 

0.162 

0.414 

0.324 

65 

0.191 

0.489 

0.383 

70 

0.237 

0.605 

0.473 

# 

75 

0.313 

0.799 

0.625 

80 

0.466 

1.191 

0.932 

85 

0.928 

2.373 

1.857 
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